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We present a reaction/diffusion model for the formation of a lower polariton condensate in a micro 
cavity containing an organic semiconducting molecular crystalline film. Our model-based upon an 
anthracene film sandwiched between two reflecting dielectric mirrors-consists of three coupled fields 
corresponding to a gas of excitons created by an external driving pulse, a reservoir of vibron states 
formed by the coupling between a ground-state vibrational model and a cavity photon, and a lower 
polariton condensate. We show that at finite temperature, the presence of the vibron reservoir can 
augment the exciton population such that a lower critical pumping threshold is required to achieve 
condensation. 



I. INTRODUCTION 

Over the past few years, there has been considerable 
excitement regarding the formation of polariton Bose 
condensates in quantum dot microcavitiesPELD This has 
spurred a parallel effort to observe similar effects in 
micr ocavity systems involving organic semiconducting 
films !il!2H2D While organic semiconducting systems share 
many of the features as their inorganic counterparts, they 
differ in that excitons tend to be more local to the individ- 
ual molecular sites, phonons play a central role, and the 
dielectric constant is considerably smaller. As a result, 
organic systems present both a challenge and opportu- 
nity for observing optical/electronic processes! 17 * 18 ^ 

In this paper we consider the formation of polariton 
condensates under non-equilibrium conditions. While 
Bose condensation in atomic gases can be understood 
as an equilibrium phenomona since one can impose the 
condition that the average number of atoms in the gas 
is conserved, polaritons form due to the strong coupling 
between a photon field and an optical transition. They 
are transient quasi-particles that can exist only so long 
as the system is continuously resupplied with photons. 
Bose condensation occurs once the density of excitons 
reaches a critical threshold such that the condensate den- 
sity grows exponentially in time. 

For the cases in consideration here, we assume that ex- 
citons and photons within the cavity couple to give rise to 
at least two incoherently coupled populations: free exci- 
tons and lower polaritons. We can also add to this picture 
quantized vibrational modes of the molecular electronic 
ground. Let us denote the vibronic states as \n v ) where 
n denotes electronic excitation and v denotes vibrational 
excitations. In the weak coupling limit, optical transi- 
tions occur predominantly from the lowest vibrational 
state of the electronic ground state to various vibrational 
states in the excited electronic state, with probability 
weighted by the Franck-Condon factors. However, in the 
limit of strong coupling to the photon field, which is ap- 



propriate for discussing polaritons, the lower electronic 
state is dressed by the photon field and vibrational states 
in the electronic ground comes into resonance with the 
states in the upper electronic state. Here we consider a 
three-state system in which the ground state, |0 o ), and 
a single quantized mode of the ground state, |0i), are 
dressed by the photon field and are strongly coupled to 
the |1 ) state in the upper electronic state. Thus, the 




FIG. 1: Dispersion curves for a system with an exciton (with 
frequency ui x ) and vibrational mode (with frequency uj v ) is 
coupled to a 2D photon cavity. The curves are labels to re- 
flect the dispersion of the lower polariton (LP) and two upper 
polariton (UPi and UP2) branches. Here the cavity off-set, 
A, is such that uj x = A + ui v and the dashed line indficates 
the bare exciton energy. 
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dressed states include the photon field in the sense that 

(0 o ,n k + l\H\0 o ,n k + l) = hu c {k) (1) 
(0i,n fc + 4ff|0i,n fc + l) = h{uj v +uj c {k)) (2) 
(l ,n k \H\l ,n k ) = Tluj x , (3) 

where u> c (k) ~ k 2 j 2m + A is the dispersion relation for 
a microcavity with cut-off frequency A = 2irc/ ' Lrj and ef- 
fective mass m = 2 r Kr\j cL. Here, r\ is the refractive index 
of the material in the cavity and L is the spacing between 
reflective mirrors. For example, a L =2300A cavity con- 
taining anthracene (77 = 1.70), one obtains h/S. = 3.11eV 
and an effective mass of m = 1.8 x 10~ 5 m e . A cavity with 
this geometry would be in resonance with the S a —> Si 
electronic transition of anthracene. This amounts to a 
"shift" in the lower energy levels by A. In other words, 
the presence of a photon in the cavity can bring any of the 
vibrational states into resonance with the exciton when 
the cavity itself is not in resonance with the excitons. 

The coupling between the exciton state and the cavity 
states is given by the Rabi frequencies 

ftfii = (l ,n fc |ir|Oo,n fc + l) 

and 

nn 2 = (lo,n k \H\Oi,n k + l) 

as given by the transition moment between the |lo) ex- 
citon state and the two vibrational states in the ground 
electronic state, |0 ) and |0i). In general, the Rabi fre- 
quency for an on-resonance excitation can be expressed 
in terms of the electric field amplitude of the cavity mode 
and the transition dipole moment viz. 

n = vfiE/n (4) 

If we invoke the Franck-Condon principle and assume 
that the transition occurs in a fixed frame of the nuclear 
motion, then 

M/i = Mio/o-w (5) 

where fo- v < 1 is the overlap integral between displaced 
harmonic oscillator states and /zio is the exciton transi- 
tion moment. Thus, knowing the Rabi frequency, f2i for 
the exciton in the cavity we can reliably estimate 

^2 = /o-i^i = /o-i(/o-oMio-E/ft)- (6) 

The overlap is determined by the Huang-Rhys factor, S 
with fo-v = v S v /v\ exp(— S). For polyacene systems 
such as anthracene, the Huang-Rhys factor for the high- 
frequency C=C modes between the v = vibrational 
state of the exciton and the v = 1 vibrational state of 
the ground electronic state is approximately 5 = 1 for 
the lowest excited state as d eterm ined by comparing the 
(0-0) and (0-1) spectral areas! 121221 This implies that fl 2 ~ 

0.36 x n x . 

Diagonalizing H in the dressed basis, one obtains po- 
lariton dispersions similar to what is shown in Fig. [l] 



Here, we have taken the exciton energy at 3.11 eV, which 
is the S — > Si transition in anthracene, included a sin- 
gle vibrational mode at 0.2 eV, and set the cavity cut-off 
at 2.9 eV so that A = uj x — u> v . The resulting polariton 
states are superpositions of a cavity photon mode with 
the molecular exciton and ground-state vibrational mode 
(which we'll refer to as a "vibron" mode). The lower po- 
lariton (LP) branch, which is predominantly "photon" in 
character undergoes avoided crossings with the two up- 
per polariton branches which are predominantly vibronic 
or excitonic in character. 



A multi-state model involving the mixing between the 
cavity and the vibronic sub levels was used recently by 
Kena-Cohen and ForreslPS^Il to describe the multiple res- 
onances observed in anthracene microcavites. In this 
case, the vibronic sublevels are those of the electronic ex- 
cited state and not the ground state, which we assume. 
At thermal equilibrium, all of the population is in the |0q) 
state and the vibrationally excited states would have es- 
sentially no population and would not ordinarily come 
into play in the optical response of the system. However, 
in a driven system and one in which the | lo) exciton state 
does carry oscillator strength to the |0„) states via the 
Franck-Condon principle, one should be able to create a 
non-thermal population inversion in these states which 
can act as a reservoir. 



In this paper, we consider the role that vibrational 
cooling may have on the formation of a polariton conden- 
sate in the incoherent pumping regime. In a recent set of 
experiments by Lidzey's group, it was apparent that vi- 
brational modes can contribute to the LP polariton popu- 
lation and may serve as a n imp ortant cooling mechanism 
for UP and free excitons! 15 * 16 ! Our model consists of a 
reaction/diffusion equation whereby a population of free 
excitons is coupled to a population of vibrational modes 
and can serve as a reservoir for the formation of a LP 
Bose condensate. Our model assumes that condensation 
occurs when fluctuations about a steady state condensate 
with zero initial amplitude begin to grow exponentially 
in time once a threshold population of excitons has been 
achieved. The model results in a series of phase-diagram 
of exciton pumping rate vs. lattice temperature and cav- 
ity cut-off. We show that when the cavity cut-off reso- 
nance from the exciton mode, thermalized vibrons can 
enhance the exciton population and lead to lower critical 
pumping rates as the lattice temperature is increased. 
This seems counterintuitive since one expects conden- 
sates to form as the temperature is lowered. However, 
as we discuss next, the symmetry breaking transition in 
a driven system comes about due to the external driving 
force rather than through internal interactions or density 
of states. 
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II. THEORETICAL REACTION/DIFFUSION 
MODEL 

We begin with the requirement that the condensate 
density obey a continuity equation with sink and source 
terms representing the temporal decay of the condensate 
with effective mass rriLP due to cavity loss (with rate 7) 
and creation of new condensate due to the presence of 
a singlet exciton reservior, S. Thus, we write using the 
hydrodynamic forrrP^l 

d t p = -V ■ (pVs)/m LP + (r -u- j)p, (7) 

where the action, s satisfies a quantum Hamilton- Jacobi 
equation of the form: 

^ + ^(v S ) 2 -£^Vp + ,p = o, (8) 

and that the current J = pVs/m. Using = 
v /pexp(is/7i) as the condensate amplitude, one obtains 
a Gross-Pitaevskii (GP) equation that properly incorpo- 
rates the source and sink terms as an optical potential 
term. 

l 9 t = -7 7 V 2 ^+ (gW 2 + l -( r ■ S--y))<j> (9) 

where 77 = h/2rriLp is the LP diffusion constant and g, 
r, and 7 are rates. If we take the exciton reservoir to be 
constant (in both time and space), then we can begin to 
see the conditions necessary for establishing a condensate 
in a driven, non-equilbrium system. The trivial solution 
is given by <f> = in which no condensate will forms at 
any time. This provides us with the trivial steady-state 
solution of 4> ss (x,t) — 0. 

However, if we introduce non-trivial fluctuations about 
the steady state, 4>(x,t) = <j) ss {x,t) + S(j)(x,t) and ask 
whether or not such amplitude fluctuations will grow or 
decay in time, one sees that if the optical potential term 
(r • S — 7) < 0, then any fluctuation in the polariton 
component will decay exponentially in time and no long- 
lived condensate will form. On the other hand, if the 
exciton density is sufficiently large that (r • S — 7) > 0, 
then the condensate amplitude will grow exponentially 
with a rate (r • S — 7). Thus we conclude there exists a 
critical exciton density, S Q whereby the polariton decay 
due to cavity leakage is exactly counterbalanced by re- 
population from a free exciton reservoir. Once a critical 
density of excitons has been achieved in a given region, 
the LP condensate population in that region experiences 
exponential growth. Thus, one expects that will occur 
once S = 7/r. Taking the exciton density to be propor- 
tional to the exciton pumping rate, p G , then one expects 
Po,crit oc 1/r. The presence of a threshold exciton pop- 
ulation for the formation of a condensate at long time 
is the hallmark of a quantum phase transition in a non- 
equilibrium system! 2 " 26 ! The principle result of the work 
presented here is that the inclusion of dressed ground 
state vibrational modes 



As discussed in the Introduction, local vibrational 
modes play very important roles in the electronic spec- 
troscopy of most molecular semiconductors and in micro- 
cavity system, may serve as a viable cooling mechanism 
for excitons. With this in mind, we adopt a continuum 
model where we write the singlet exciton density, S(x, t), 
and vibron density, v{x,t), as a solutions to the following 
reaction/diffusion equations which, in turn, are coupled 
to the GP equation for the condensate amplitude, 

d t S = D S W 2 S -r ■ S\c/)(x,t)\ 2 +p(x,t) 

- j s S+(k lV (x,t)~kS) (10) 
d t v = D v V 2 v - k lV {x,t) + kS (II) 

where r>S|(/)| corresponds loss of exciton density in a 
given spatial region due to the formation of the conden- 
sate, 75 is the singlet exciton radiative rate, and p(x, t) 
is the rate that excitons are pumped into the system by 
some external source. We imagine a experimental situ- 
ation where excitons in a thin-film within a microcavity 
are created in a given region illuminated by a steady laser 
spot, the excitons diffuse incoherently through through 
the sample as UP polaritons and are either rapidly ther- 
malized with UP vibrons or decay through radiative pro- 
cess. Taking A = na as the area illuminated by the 
pumping pulse we can relate the pumping rate to an ir- 
radiance via E e = p a uj x /A ~ 0.32p o W/m 2 taking a — 
500/im 2 and p is given in ns _1 . Thus a p Q = 4000ns -1 
would correspond to an irradiance of 1.3 kW/m 2 which is 
comparable to the average solar irradiance on the Earth's 
surface. 

Putting aside any decay channels for the moment, let 
us assume that the exciton and vibron populations are 
rapidly thermalized so that we can write 

fci = e~ pt k 

where S = uj x — (lu v + A) is the energy gap between the 
exciton and vibrons. This provides a useful means for 
incorporating the geometry of the cavity and the tem- 
perature of the lattice itself into our equations of motion. 
We also have assumed that the energy difference between 
the free exciton reservoir and LP polariton condensate is 
much greater than fc^T so that any exciton population 
removed by condensate formation does not escape to re- 
populate the free exciton reservoir. A final assumption 
we make is that the diffusion constant appearing in Eq. |ll| 
can be effectively set to zero since we can assume that 
molecular vibrational excitations are local to the molec- 
ular sites in the system. 



III. NUMERICAL RESULTS 

We next consider the numerical integration of the equa- 
tions of motion given above. For this, we use the method 
of lines approach as imp lemented in the NDSolveQ rou- 
tine in Mathematical For numerical purposes, we use 
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the exciton decay rate, 75 and the singlet exciton dif- 
fusion constant to set the time and length-scales of our 
model. The exciton decay rate, 75 includes all radia- 
tive and non-radiative processes that diminish the exci- 
ton population that do not contribute to the formation 
of either lower polaritons or |0i,n + 1) vibrons. Since the 
singlet exciton radiative lifetime is at least 1 ns, setting 
75 = Ins -1 is a reasonable estimate. From transient 
singlet diffusion constant for e xcitons in anthracene is 
D s = 0.13^m 2 /ns at 1 OK .1221221 This produces a length- 
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excitons are 
= IOOtt^ 2 , so 



scale I = \jDsllS = 1.14/xm. From Eq. 
pumped into the system in an area A oc ira 
that poj 'A gives the number of excitons that are created 
in the system per unit time per unit area. 



Below, we consider two possible regimes. One in which 
lifetime of the exciton gas is short relative to that of the 
polariton condensate and the other in which the exciton 
gas is long-lived. In each case, we "seed" the condensate 
with a small amplitude about the center of our grid and 
determine whether or not this "seed" fluctuation grows 
or decays as the temperature and pumping strength are 



varied. 



A. Slow interconversion regime 

Let us first consider the case where the conversion from 
the exciton gas to the LP condensate is slow compared 
to the exciton to vibron conversion. We shall refer to 
this as a slow conversion regime and take r « k in 
our equations of motion. We would expect this limit to 
be valid when the energy difference between the vibron 
state and the free exciton is on the order of kT and the 
Rabi splitting between the exciton and cavity modes to 
be small. Taking the pumping rate to be a gaussian that 
grows exponentially to some maximum value 



p(x,t)=p e-* 2 ^(l-e- t / T ) 



(12) 



In Fig|2]we show the results of a representative numer- 
ical simulation using a pumping rate of p — 4, 000ns -1 , 
which corresponds to an irradiance of s» 3.1kW/m 2 , 
which is above the threshold for forming the LP conden- 
sate. As point of reference, solar irradiance of the earth 
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is on the order of 1 kW/m 2 . The three panels in Fig[2ja- 
c) are the time-dependent densities for the free exciton 
gas, the vibron reservoir, and the condensate at a lattice 
temperature of 10K. As seen in Fig. [2^,, the exciton gas 
forms and begins populating the vibron reservoir. As the 
population in the vibron reservoir increases, it augments 
the exciton density such that a critical population can be 
established for forming a stable condensate. 

In Figs [2]i-f we show the formation of the condensate 
when the lattice temperature is considerably warmer at 
200K. Here the vibron reservoir rapidly augments the 
exciton population and the condensate forms at a much 
earlier time than in the previous example at 10K. Also, 
one notes that the condensation is far more abrupt and 
produces a rapidly evolving condensate wave packet. The 
system does, however, reach a steady state at long time. 

Once the condensate forms, it depletes the exciton den- 
sity and at long times, the system (exciton gas, vibron 
reservoir, and condensate) approaches a steady state. 
This is evidenced in Figj3[a-c) were we show the con- 
densate density at x — (center of the grid) versus time 
for various cavities at different temperature. In each case, 
the exciton pumping rate is above the critical threshold. 

In Fig. [3^, and b, the vibron state (at A + oj v ) is lower 
in energy than the exciton energy. In both cases, the con- 
densate forms about 1 ns after the pumping is initialized. 
In all cases, the initial rapid buildup of polariton density 
drops almost immediately due to the ballistic spread of 
the polariton wave packet away from the pumping region. 
After one or two "bounces" the polarition density relaxes 
to a steady state population. 

In Fig. [3^,, the cavity off-set at A = 2.85eV is such that 
the vibron level is considerably lower than the exciton 
level. In this case, the initial formation time of the po- 
lariton condensate shows a very weak dependency on the 
temperature, but the steady-state population is clearly 
dependent upon the temperature with warmer systems 
leading to a greater steady state polariton density at long 
times. In the intermediate case at A = 2.88eV, the po- 
lariton steady state is reached much sooner. While the 
rate at which the polariton population approaches the 
steady state does depend upon the temperature, the fi- 
nal steady-state population shows very little dependency. 

Lastly in Fig. [3j; we consider a case in which the vibron 
level is well above the exciton level and plays no signifi- 
cant role in the formation of the polariton condensate at 
any temperature. Here, we can again see a series of beats 
as the polariton condensate forms, decays, and reforms 
until a steady-state population is achieved. These beats 
are not due to quantum coherences between the exciton 
and polariton since we are treating the exciton gas as 
a classical reservoir. The beats are are due to popula- 
tion depletion and replenishment from both the vibron 
reservoir and the pumping field. 

In examining the final steady state populations, we 
note that steady-state condensate population is enhanced 
by at least a factor of 3 to 5 by the presence of the vi- 
bron reservoir. This implies that cavities which can take 



advantage of the lower-lying vibron levels can produce a 
higher density condensate which translates into a more 
intense steady-state optical signal. 

In Fig. [4] we show how the condensation delay times de- 
pend upon both the lattice temperature and the exciton 
pumping rate for a cavity with A = 2.88eV. The model 
does give a critical pumping rate of p a = 3800ns -1 , cor- 
responding to an irradiance of 1.2 kW/m 2 , below which 
long-lived polariton condensates will not form at any 
temperature. Close to this threshold, the delay time is 
very sensitive to the lattice temperature to the extent 
that for a very cold lattice the condensation delay time 
can be rs 6x longer that of a warmer lattice. It is impor- 
tant to point that an irradiance of 1.2 kW/m 2 is essen- 
tially the intensity of bright sunlight which implies that 
the critical pumping rates given by our model are not un- 
reasonably intense and can easily be achieved. However, 
this is the irradiance within the cavity itself and we have 
not taken into account the fact that the pumping laser 
needs to penetrate through the DBR to the cavity. 

On first thought it seems counterintuitive that a 
warmer system would have an easier time forming a sta- 
ble polariton condensate since we expect condensation 
to occur as the temperature is lowered past some critical 
temperature. However, the symmetry breaking mecha- 
nism in this case is the density of excitons in a given 
region which are introduced into the system via the ex- 
ternal pumping source. By populating the vibron mode, 
which does not contribute directly to the lower polariton 
population, and allowing this mode to equilibrate and 
populate the exciton mode, it effectively augments the 
local exciton density such the the critical density can be 
reached at lower exciton pumping rates. Increasing the 
temperature, causes as shift in the steady state popula- 
tions of vibrons and excitons towards more population in 
the higher-energy species. 

B. Rapid interconversion regime 

Central to our theory is the rate that free excitons are 
converted to lower polaritons. To estimate this rate, let 
us assume that the "golden rule" is valid, such that the 
maximum value of the rate is given by 

r max = ^M 2 p{k ex = k cav ) (13) 

where p(k ex = k cav ) is the density of cavity states eval- 
uated at the avoided crossing. Since the cavity con- 
tains a two-dimensional gas of photons with dispersion 
Ek = A + fink' 1 , the density of states is constant with 

P(k ex = k cav ) = X- (14) 

nn 

Using the numerical values for our model cavity and set- 
ting = 0.05eV, one obtains r max = 1.9 x lO'ns . 

We can follow a similar line of reasoning for the cou- 
pling between the exciton and vibron reservoirs. The vi- 
bron reservoir in our model corresponds to an ensemble 
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FIG. 3: Condensate density |</>(t)| 2 at a; = for different cavities for p above the critical pumping threshold. 



of UP polaritons formed by mixing a cavity mode with a 
ground- electronic state vibrational mode of the molecules 
in the cavity. If we assume that the golden rule holds, 
then we can estimate the exciton to vibron conversion 



rate by multiplying Eq. 13 a Franck-Condon factor, i.e. 



x rv 



Based upon arguments above, the ex- 



citon to vibron conversion rate is at least an order of 
magnitude slower than the exciton to LP conversion. 

For the LP decay, let us assume that the cavity has a 
quality factor, Q — 6000, that corresponds to the number 
of times a photon will traverse the cavity before escaping. 
Taking the cavity width to be w 250 nm, one obtains an 
LP decay rate of 7 = 120ns -1 . 

Since the conversion rates in this regime are far greater 
than in the slow conversion regime we just examined, 
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FIG. 4: Delay of condensate formation in a "standard" an- 
thracene cavity model. (A = 2.88eV, U) x = 3.1eV, and 
u v = 0.2eV). 
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FIG. 5: Condensate amplitude at x = for a non-resonant 
cavity in the "rapid-conversion" model following a 30fs pulse 
with peak intensity of 60 ns -1 . Each curve is labeled by the 
lattice temperature (in K). 



the critical pumping strengths will be considerably lower. 
However, the life times and dynamics will be consider- 
ably faster. In thinking about possible experimental sit- 
uations, we consider the effect of a single 30 fs gaussian 
shaped pulse that drives the creation of an exciton gas. 
The density of excitons created will depend upon the in- 
tensity (and duration) of the pulse. Prior to the arrival 
of the pulse, we seed the condensate with a very small 
amplitude fluctuation at x — which will have decayed 
somewhat by the time the pulse arrives. 

In Fig. [5] we show how the condensate population at 
the center of the simulation cell varies in time for a non- 
resonant cavity system with A = 2.88eU and at a hxed 
pumping strength of p = 60ns -1 . At t = 0/s, a conden- 
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FIG. 6: Phase diagrams for non-resonant cavities in the 
"rapid-conversion" model. 



sate fluctuation is introduced and at t = 60fs the 30 fs 
excitation pulse reaches its maximum intensity and cre- 
ates a population of excitons, which go on to produce 
more condensate as well as vibrons. After 100 fs, the 
system evolves without further pumping and the popu- 
lation in the condensate will either decay due to cavity 
loss, or continue to grow as excitons are converted to LP 
condensate. The various curves shown here correspond 
to different lattice temperatures. For the case at hand, 
when the lattice temperature is above 45K, the conden- 
sate population continues to grow following the initial 
pulse. This indicates that sufficient exciton density is 
present to spawn the formation of a long-lived polariton 
condensate. Below this temperature, the LP population 
simply decays. 

Fig.|6]shows a phase diagram for the off-resonant cavity 
in terms of the threshold temperature needed to produce 
a long-lived condensate at a given pumping strength. As 
temperature is decreased, the role of the vibronic reser- 
voir is diminished and higher pumping strengths are re- 
quired to from a long-lived condensate. 

In Fig. [7] we considered what happens if there is a spa- 
cial off-set between the pumping pulse and the initial 
condensate "seed" fluctuation. We initialized the con- 
densate with a small gaussian centered off-center at 50ttm 
and created exciton density centered about x = at 60 
fs later. In this case, the condensate appears to migrate 
before decaying and there are a number of "beats" in 
the density itself. The beats occur even when the seed 
fluctuation is located at x = and are due to cyclic 
growth/decay kinetics between the condensate and ex- 
citon densities, much like what is observed in preda- 
tor/ prey and disease models rather than to coherent mo- 
tion of the condensate itself. At long times following 
the pulse, the condensate eventually decays due to cav- 
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FIG. 7: Time evolution of condensate in the "rapid conver- 
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ity loss. It should be pointed out, however, that while the 
requisite pumping intensity in this model is lower than in 
the "slow conversion" case, the integrated density of the 
condensate wave function is considerably less. However, 
this density can be used to further seed the production 
of additional condensate if subsequent pulses are intro- 
duced. 
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the Bogoliubov quasi-particle spectra for the equilibrium Bose 
condensate. Unless otherwise specified, we have set, g = 1, 
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1 and rS = 3. 



IV. STABILITY ANALYSIS 

We can introduce a time-dependent fluctuation about 
the steady-state ^> = {(f> , ip*,S , v a } by writing the so- 
lution to Eqs. 9 and 10 as *(0) = * + 5V(t) 

Let us consider the time evolution of Eqs. 9 and 10 by 
writing $>(f) = {</>(*), if>*(t), S(t), v(t)} and introducing a 
fluctuation about its steady state such that ^f(t) = f + 
S^f(t). Introducing this into Eq. 9. and 10, keeping only 



terms linear in S^S, and transforming from real to k-space, 
produces a set of coupled equations for the fluctuations: 



/ <50k \ 


t 




= M(k) 


V Sv k J 


V 



where the stability matrix is given by 



M(k) 



7 ? k 2 + ( 5 |^| 2 + i(r5 - 7 )) 

-i(f>*S 




-?yk 2 



g(M 2 

(5l^o| 2 -|(r5 - 7 )) 
-i<j) S 




%r<f> 










+ k + r\(j) \ 2 ) 


ik' 


ik 


-ik 



(16) 



In writing this, we assume the steady-state solutions are 
homogeneous and that wave-vector, k, remains a con- 
stant of the motion. We have also used fj, = g\<p \ 2 as the 
chemical potential of the condensate. ( Note that the 
k and k' appearing in the matrix elements coupling the 
reservoirs are rate constants and not wave- vectors.) The 
stability matrix, A4, given above is similar to the equa- 
tions of motion reported by Wouters and CarusottcE^ 
and by Byrnes et alW^ In the absence of the reservoir, 
the eigenvalues of the upper 2x2 block of M(k) are 
the Bogoliubov modes corresponding to excitations of the 
condensate. These modes vanish at k = correspond- 
ing to a Goldstone brach that can be understood as a 
slow rotation of the condensate across the sample^ and 
the linear dispersion about k = is characteristic of a 
superfluid state. 

Let us consider the stability of the polari- 
ton/exciton/vibron equations in the steady-state 
limit by first assuming that the exciton and vibron 
reservoirs are equilibrated rapidly. This allows us to 
write an effective exciton decay constant, 7$, that will 
depend both upon the lattice temperature and the 
cavity offset. Including the coupling to the reservoir, 
two regimes can be identified. First, in the limit where 
7^ 3> 7 in which the reservoir relaxation is very rapid 
compared to the polariton relaxation rate, the disper- 
sion around k = is flat corresponding to a diffusive 
Goldstone mode which is in stark contrast to the linear 
dispersion of the sound-mode in an equilibrium Bose 
condensate. For a Bose condensate, the gap at k = 
necessarily vanishes and the dispersion about k = is 
linear corresponding to the formation of a superfluid 
with the slope proportional to the sound velocity. In 
the non-equilibrium case, the sound velocity vanishes as 
well as seen in Fig. [8^, which corresponds to our "slow" 
regime in which the exciton gas decays rapidly compared 
to the polariton condensate. 

On the other hand, if 7 ^> j' s a gap opens at k = 



and the polaritons again take on an effective mass. The 
gap and the flattening of the spectrum is akin to the 
blue-shift observed by Kasprzak for a polariton conden- 
sate driven well above the condensation thresholdP^Ult 
was speculated that the flattening could be due to either 
polariton localization in real-space or to pump and decay 
processes. Here, as in Refs l24l36l the flattening and gap 
is entirely due to competing decay processes between the 
exciton reservoir and the polariton gas. 



V. DISCUSSION 

In this paper we have focused upon the role that 
molecular vibrational excitations may play in the forma- 
tion of a LP condensate in a microcavity containing a 
thin-film organic semiconductor, such as anthracene. In 
our model, molecular vibrational excitations coupled to 
the photon field serve as a secondary thermal reservoir 
for free excitons. Thermal fluctuations from the vibron 
reservoir can augment the exciton density such that at 
finite lattice temperature polariton condensation can be 
achieved at lower exciton pumping rates. We assumed 
that the vibron reservoir involved only a single vibra- 
tional model per molecule. In a realistic system, there 
will be multiple vibrational modes as evidenced by the 
various Franck-Condon peaks in the emission spectrum 
of many organic semiconductors. Assuming the mecha- 
nism presented here holds true in the multi-mode case, 
these Franck-Condon modes should provide a reservoir 
for the exciton gas and greatly facilitate the formation of 
a stable LP condensate. 

We explored two limits of the model. In the "slow 
conversion" model, we assumed that the rate of conver- 
sion between continuously driven and replenished gas of 
free excitons and a vibronic reservoir level was compara- 
ble to the rate of conversion between the exciton and a 
lower polariton condensate. In this limit, the tempera- 
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ture of the molecular crystal lattice plays an important 
role in providing a sufficient exciton density to spawn the 
growth of fluctuations about the condensate vacuum. At 
long times, this model produces steady state solutions 
with populations in both condensate and non-condensate 
channels. In the "rapid conversion" limit, we assumed 
that the conversion of the exciton gas to LP condensate 
could be estimated from the Rabi splitting between the 
upper and lower polariton branches and the density of 
states of photons in a two-dimensional cavity. In this 
model, the exciton gas is pumped impulsively via a single 
30 fs pulse. Here, too, the growth or decay of conden- 
sate fluctuations hinges upon the pumping being above 
a threshold intensity and shows a pronounced sensitivity 
to the temperature of the lattice when the cavity is off- 
resonance with the exciton. In both models, it is highly 
encouraging that the threshold intensities are well within 
the margins of current experimental capabilities. 



Lastly, along a similar vein, one can model an electri- 
cally pumped system in which electrons and holes are in- 
jected onto a microcavity device and recombine to form 
singlet and triplet excitons. We suspect that both the 
triplet excitons and bound charge-transfer contact pairs 
may serve a similar role as the vibronic reservoir in the 
present model. We are currently exploring this possibil- 
ity. 
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